UCI-TR- 201 1-03, UH511-1157-2011 



Isospin- Violating Dark Matter 



O 
(N 
^— » 
O 

O 



Oh. 
> 

m 

O 



X 



Jonathan L. Feng/ Jason Kumar, ^ Danny Marfatia,^'^ and David Sanford^ 

^Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA 
"^Department of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822, USA 
^Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA 
"^Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

Searches for dark matter scattering off nuclei are typicahy compared assuming that the dark 
matter's spin-independent couphngs are identical for protons and neutrons. This assumption is 
neither innocuous nor well motivated. We consider isospin- violating dark matter (IVDM) with 
one extra parameter, the ratio of neutron to proton couplings, and include the isotope distribution 
for each detector. For a single choice of the coupling ratio, the DAM A and CoGeNT signals are 
consistent with each other and with current XENON constraints, and they unambiguously predict 
near future signals at XENON and CRESST. We provide a quark-level realization of IVDM as 
WIMPless dark matter that is consistent with all collider and low-energy bounds. 

PACS numbers: 95.35.+d, 12.60.Jv 



Introduction. Dark matter makes up five-sixths 
of the matter in the Universe, but all current evi- 
dence for dark matter is through its gravitational ef- 
fects. The detection of dark matter scattering through 
non-gravitational interactions would be a large step to- 
ward identifying dark matter, and there are many ex- 
periments searching for such events. The excitement 
around this approach has been heightened recently by 
data from the DAMA [D] and CoGeNT [H experiments, 
which are consistent with scattering by a dark matter 
particle with mass mx ~ 10 GeV and spin-independent 
(SI) X-nucleon scattering cross sections (Jn ^ 2x 10~^ pb 
and 5 x 10 ~^ pb, respectively. This excitement is, how- 
ever, tempered by null results from XENON [l, Q and 
CDMS jsl, Q , leaving a confusing picture that has moti- 
vated much theoretical and experimental work. 

The comparison of dark matter experimental results 
is subject to an array of assumptions and uncertainties 
from particle physics, nuclear physics, and astrophysics. 
In this study, we focus on a particularly simple and com- 
mon particle physics assumption, that of flavor isospin 
invar iance. Dark matter detectors have various nuclear 
compositions. To derive implications for (Jat, experi- 
ments almost universally assume that dark matter cou- 
ples identically to protons and neutrons. This assump- 
tion is not well motivated. For example, Dirac neutrinos 
and sneutrinos have isospin-violating couplings, and, in 
fact, even neutralino couplings are generically isospin- 
violating, although typically at insignificant levels (see, 
e.g., Refs. 0, i). In any case, as these and other con- 
ventional dark matter candidates do not easily explain 
the DAMA and CoGeNT signals, it is reasonable to con- 
sider more general frameworks. Here we consider flavor 
isospin-violating dark matter (IVDM) with one extra pa- 
rameter, the ratio of neutron to proton couplings fn/fp- 

IVDM has been considered previously in general anal- 
yses [sl, and also recently in studies of various interpre- 
tations of the CoGeNT results [iQ|,[il|. We focus solely 



on IVDM and consider for the first time the distribution 
of isotopes present in each detector. Previous work has 
neglected this distribution, which implies that dark mat- 
ter may be completely decoupled from any given detector 
for a particular value of fn/fp- However, this is not true 
if there is more than one isotope present, as is the case 
in many detectors, and the viability and implications of 
IVDM cannot be established without considering the iso- 
tope distribution. As we will see, including the isotope 
distribution has remarkable consequences. For a single 
choice of /n//p, the DAMA and CoGeNT signals are 
consistent with each other and with current XENON con- 
straints. At the same time, the isotope distribution im- 
plies that XENON cannot be completely decoupled, and 
the IVDM scenario unambiguously predicts near future 
signals at XENON and other detectors, such as CRESST 
and COUPP. We identify and discuss slight inconsisten- 
cies with other data, and present a general analysis of 
when experiments may be reconciled by isospin violation. 
Finally, we provide a quark-level realization of IVDM as 
WIMPless dark matter p^, ll3| that is consistent with all 
collider and low-energy bounds. 

Cross sections for IVDM. We focus on the SI scat- 
tering of an IVDM particle X off a nucleus A with Z 
protons and A — Z neutrons. The event rate is 



R NtUx 



(fvf{v)v 



da 



(1) 



where Nt is the number of target nuclei, nx is the lo- 
cal number density of dark matter particles, and the 
limits of the recoil energy integral are determined 
by experimental considerations. The IVDM particle's 
velocity v varies from Vmin = \/'^aEr/2ii\, where 
fiA = "ruATTix / {tua + TTix)^ to '^max, a fuuctiou of the 
halo escape velocity, and f{v) is the distribution of X 
velocities relative to the detector. The differential cross 
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section is da/dEn = aATnA/{'2v'^fi\), with 

cta = ^[UZF^AEr) + MA-Z)F2{En)f , (2) 

where fp^n are the coupHngs to protons and neutrons, 
normahzed by the choice of mass scale M*, and F^^{Eb) 
are the proton and neutron form factors for nucleus A. 

FI{Er) and F'1{Er) are not identical. FI{Er) is 
what has typically been measured, but F^{Eb) may also 
be probed, for example, through neutrino and electron 
parity- violating scattering off nuclei Q. However, since 
the isospin violation from this effect is small compared to 
the potentially large effects of varying /n//p, we will set 
both form factors equal to Fa{Er). With this approxi- 
mation, the event rate simplifies to R = ctaIa^ where 

^A = §-JUZ^U{A-Z)f (3) 
Ia = NtUx I dEn f^^^d\ f{v)^Fl{ER) , (4) 

and aA is the zero-momentum-transfer SI cross section 
from particle physics, and Ia depends on experimental, 
astrophysical, and nuclear physics inputs. If fn = /p, 
we recover the well-known relation R (x A^. For IVDM, 
however, the scattering amplitudes for protons and neu- 
trons may interfere destructively, with complete destruc- 
tive interference for fn/fp = —Z/{A — Z). 

We assume that each detector either has only one el- 
ement, or that the recoil spectrum allows one to distin- 
guish one element as the dominant scatterer. But it is 
crucial to include the possibility of multiple isotopes. The 
event rate is then R = ^iTjiCJAjAi^ where the sum is 
over isotopes Ai with fractional number abundance rji. 

IVDM and current data. It will be convenient 
to define two nucleon cross sections. The first is <jp = 
/ip/^/M^, the X -proton cross section. In terms of a^, 

i 

The second is a^, the typically- derived X -nucleon cross 
section from scattering off nuclei with atomic number 
Z, assuming isospin conservation and the isotope abun- 
dances found in nature. With the simplification that the 
lAi vary only mildly for different z, we find 

If one isotope dominates, the well-known result, Fz = 
[Z/A -h (1 - is obtained. 

In Fig. [T]we show regions in the (mx,cr^) plane and 
the {mx^cTp) plane for fn/fp = —0.7 that are favored and 
excluded by current bounds. These include the DAM A 
3<j favored region [Til, [l6| , assuming no channeling [VJl 
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FIG. 1. Favored regions and exclusion contours in the 
{mx,crfi) plane (top), and in the {mx,crp) plane for IVDM 
with fn/fp — —0.7 (bottom). 



and that the signal arises entirely from Na scattering; the 
CoGeNT 90% CL favored region (2|; 90% CL exclusion 
contours from XENONIOO [sj and XENONIO and 
90% CL bounds from GDMS Ge and Si [5, 6]. The isotope 
abundances are given in Tables |T] and |TT1 

There are controversies regarding the exclusion con- 
tours for xenon-based detectors at low mass [18]. The 
energy dependence of the scintillation efficiency at low 
energies is uncertain, and there are questions about the 
assumption of Poisson fluctuations in the expected pho- 
toelectron count for light dark matter. We have also not 
accounted for uncertainties in the associated quenching 
factors for Na, Ge and Si [19]. These issues can enlarge 
some of the signal regions or alter some of the exclusion 
curves of Fig. [TJ We have also not adjusted the favored 
regions and bounds to account for differences in the dark 
matter velocity distributions adopted by the various anal- 
yses, which would slightly shift the contours. 

Remarkably, for -0.72 < fn/fp < -0.66, the DAMA- 
and CoGeNT-favored regions overlap and the sensitivity 
of XENON is sufficiently reduced to be consistent with 
these signals, since this choice of fn/fp leads to nearly 
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TABLE I. i?max[^i, Z2], where the Zi (Z2) elements are hsted 
in rows (columns). Elements with one significant isotope have 
their (Z, A) listed; those with more than one are denoted by 
asterisks and listed in Table Ull 



Element 


Xe 


Ge 


Si 


Ca 


W 


Ne 


C 


Xe(54,*) 


1.00 


8.79 


149.55 


138.21 


10.91 


34.31 


387.66 


Ge (32,*) 


22.43 


1.00 


68.35 


63.14 


130.45 


15.53 


176.47 


Si (14,*) 


172.27 


30.77 


1.00 


1.06 


757.44 


1.06 


2.67 


Ca(20,*) 


173.60 


31.53 


1.17 


1.00 


782.49 


1.10 


2.81 


W(74,*) 


2.98 


13.88 


177.46 


166.15 


1.00 


41.64 


466.75 


Ne(10,*) 


163.65 


28.91 


4.39 


4.09 


726.09 


1.00 


11.52 


C(6,*) 


176.35 


32.13 


1.07 


1.02 


789.59 


1.12 


1.00 


1(53,127) 


1.94 


5.51 


127.04 


118.35 


20.68 


28.92 


326.95 


Cs (55,133) 


1.16 


7.15 


139.65 


127.61 


12.32 


31.88 


355.27 


0(8,16) 


178.49 


32.13 


1.08 


1.03 


789.90 


1.13 


1.01 


Na (11,23) 


101.68 


13.77 


8.45 


8.33 


481.03 


2.27 


22.68 


Ar (18,40) 


35.38 


1.87 


45.39 


42.56 


190.26 


10.32 


119.80 


F(9,19) 


89.39 


10.88 


12.44 


11.90 


425.93 


3.05 


33.47 



TABLE IL Ai for isotopes and their fractional number abun- 
dances rji in percent for all isotopes with r]i > 1%. 

Xe Ge Si Ca W Ne C 



182(27) 20(91) 12(99) 
183(14) 22 (9.3) 13(1.1) 
184 (31) 
186 (28) 



128(1.9) 70(21) 28(92) 40(97) 

129(26) 72(28) 29(4.7) 44(2.1) 

130(4.1) 73(7.7) 30(3.1) 

131 (21) 74(36) 

132(27) 76(7.4) 

134(10) 

136 (8.9) 



complete destructive interference for the proton/neutron 
content of xenon isotopes. The possibility of IVDM 
therefore brings much of the world's data into agreement 
and leads to a very different picture than that implied 
by studies assuming isospin conservation. The CDMS 
Ge constraint marginally excludes the overlapping re- 
gion, and since CoGeNT utilizes Ge, the tension between 
CoGeNT and CDMS Ge cannot be alleviated by isospin 
violation. However, it is possible that an improved under- 
standing of CoGeNT backgrounds and the energy scale 
calibration of the CDMS Ge detectors at low energy may 
resolve the disagreement (6|, llH, |2o| . 

Predictions. Further tests of the IVDM hypothesis 
may come from other detectors. If two experiments re- 
port signals suggesting the same mx , their results imply 
an experimental measurement of 



R[Zi,Z2]=a^'/a 



N 



(7) 



R[Zi^Z2] = Fz^/Fz^ is then a quadratic equation in 
fn/fp^ the solution of which enables unambiguous signal 
predictions for other detectors. 

As timely examples, consider two current experiments. 
Preliminary results from CRESST may indicate a signal 



from scattering off oxygen [21^. fn/ fp ^ —0.7 implies 
Fz=8 ^ 44. The IVDM explanation of DAMA and Co- 
GeNT therefore predicts that CRESST will see a signal 
consistent with mx ^ 10 GeV and cr^^^ ^ 8.5 a^^^^. 
Such a cross section may in fact be consistent with 
CRESST data [ii,[22|. COUPP is a CF3I detector; its 
sensitivity to low-mass dark matter arises from C and 
F scattering. For fn/fp ^ —0.7, mx ^ 10 GeV, we 



find 



^Z=6 



4.2 erf COUPP 



would be expected to report a normalized cross section 
between these values, with the value depending on the 
relative detection power of the C and F targets. 

Relative detection prospects. Although XENON 
excludes CoGeNT and DAMA signals assuming isospin 
conservation, this is not the case for IVDM. One might 
then ask: given any signal at a detector with atomic 
number Zi, what sensitivity is required for a detector 
with atomic number Z2 to either corroborate or disfavor 
this signal, allowing for isospin violation? Maximizing 
R[Zi, Z2] with respect to fn/fp determines the factor by 
which the Z2 detector must exclude the Zi signal as- 
suming isospin conservation, such that the Zi signal is 
excluded even allowing for isospin violation. Similarly, 
maximizing R[Z2, Zi] determines the factor by which the 
Z2 detector may come up short in probing an isospin- 
conserving origin for the Zi signal, while still having the 
potential to find evidence for an isospin-violating origin. 

In Table [H we present i?max[^i, ^2], the maximal value 
of R[Zi^Z2] for all possible values of fn/fp^ for many 
materials that are commonly used in dark matter de- 
tectors. The isotope composition of elements plays an 
important role in determining Rmax[Zi^ Z2]. If the ele- 
ment Z2 is composed entirely of one isotope, then it is 
always possible to choose fn/fp so that (JZ2 = and thus 
^max[^i 7^ ^2, ^2] = oo; thcsc columus have been omit- 
ted from Table [D However, if there is more than one 
significant isotope, it is impossible to achieve exact de- 
structive interference for all isotopes simultaneously, and 
so i^max[^i 7^ ^2,^2] is finite. In particular, although 
isospin violation can weaken the bounds achieved by Xe 
and Ge detectors, we see in Table [J that these bounds can 
be weakened by at most two orders of magnitude. Up- 
coming XENON results may therefore exclude DAMA 
and CoGeNT, even for IVDM; XENON bounds already 
eliminate some of the DAMA/CoGeNT overlap region 
(Fig.[T]), and will probe the entire region if XENON sen- 
sitivities are improved by an order of magnitude. 

Isospin violation in WIMPless models. So far 

we have worked at the nucleon level. We now provide a 
quark-level theory of dark matter that generically real- 
izes isospin violation. In supersymmetric WIMPless dark 
matter models 12|,[l3, dark matter particles X freeze out 
in a hidden sector with the correct relic density and inter- 
act with the standard model through connector particles 
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Y. We consider the superpotential 
W = Y.{\lXY,,ql + XiXY^^u}, + XiiXYdM , (8) 

i 

where X is a real scalar dark matter particle, QL^UR^dn 
are standard model quarks, i labels generations, and the 
connectors Yqj^^uR,dR are 4*^ generation mirror quarks. 
Assuming real Yukawa couplings and ttiy = ^Y^^d ^ 
Trix^rnq, the connector particles induce the SI operators 



{9) 



leading to the scattering cross section of Eq. (|3]) with 



-XlXl^B^f)/{^/7^mxmY)• The 
B^l^ are integrated nuclear form factors, including = 



BJ«6,i?: = B^«4 [7]. 



The amount of isospin violation in dark matter- nucleus 
interactions is solely determined by the Yukawa flavor 
structure. There are many possibilities; WIMPless mod- 
els may explain the DAMA signal with couplings to either 



the [l2| or 3^^ [l3|, |23| generation. Here we assume 



only 1®* generation quark couplings, automatically sat- 
isfying flavor constraints. Assuming mx = 10 GeV and 
my = 400 GeV, consistent with all collider and precision 
electroweak bounds, the region of the (AjA^, AjA^) plane 
that explains DAMA and CoGeNT is 



A,, 



-1.08 A^, 



0.013 < XlX^ < 0.024 . (10) 



IVDM is clearly generic in this microscopic model of dark 
matter interactions and may simultaneously reconcile the 
DAMA and CoGeNT signals and XENON bounds. 

The IVDM reconciliation of DAMA, CoGeNT, and 
XENON relies on cancellations between p and n cou- 
plings, and so requires larger couplings than in the 
isospin-preserving case to maintain the desired DAMA 
and CoGeNT signals. Such models may potentially vi- 
olate collider constraints, which are not subject to can- 
cellations. This WIMPless model provides a quark- level 
framework in which one may investigate this question. 

The most stringent model-independent constraints are 
from Tevatron searches for pp XX + jet 0, HH. Us- 
ing MadGraph/MadEvent 4.4.32 [26], one can compute 
the monojet cross section (requiring jet Et > 80 GeV) 
induced by the operator of Eq. (|9]). The resulting 2a 
bounds from Tevatron data are roughly AJAJ^ ^ ^ 1, two 
orders of magnitude too weak to probe the DAMA and 
CoGeNT favored couplings described in Eq. ([TOj). 

Conclusions. Results for spin-independent dark 
matter interactions typically assume identical couplings 
to protons and neutrons. Isospin violation is generic, 
however, and we have shown that IVDM with fn/ fp ^ 
—0.7 may explain both DAMA and CoGeNT, consis- 
tent with XENONlO/100 bounds. This scenario is only 
marginally excluded by CDMS Ge constraints, unam- 
biguously predicts a signal at CRESST, and may even be 
tested by XENON, given its several significant isotopes. 



as discussed above; near future data will shed light on this 
picture. More generally, we have explored the extent to 
which dropping the fp = fn assumption may reconcile re- 
sults from various detectors, stressing the important role 
played by the distribution of isotopes. Finally, we have 
shown that IVDM is easily realized in a quark- level model 
consistent with all low-energy and collider observables. 
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Note added. After the completion of this work, an 
annual modulation signal from CoGeNT and a new con- 
straint from SIMPLE have been reported. These results 
and some of the following discussion may be found in 
Refs. \2^,M- 
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